Preprint typeset in JHEP style - HYPER VERSION 



Holographic Hadrons in a Confining Finite 
Density Medium 



Yunseok Seo* Jonathan P. Shock; Sang-Jin Sin and Dimitrios Zoakos § 

* Center for Quantum Spacetime, Sogang University, Seoul 121-742, Korea 

' § Departamento de Fisica de Particulas, Universidade de Santiago de 
Compostela and Instituto Galego de Fisica de Altas Enerxias (IGFAE) 
E-15782 Santiago de Compostela, Spain 

^ Department of Physics, Hanyang University, Seoul 133-191, Korea 



Abstract: We study a sector of the hadron spectrum in the presence of finite baryon 
density. We use a non-supersymmetric gravity dual to a confining guage theory which 
exhibits a running dilaton. The interaction of mesons with the finite density medium 
is encoded in the dual theory by a force balancing between flavor D7-branes and a 
baryon vertex provided by a wrapped D5-brane. When the current quark mass m q is 
sufficiently large, the meson mass reduces, exhibiting an interesting spectral flow as 
we increase the baryon density while it has a more complicated behaviour for very 
small m q . 

KEYWORDS: AdS/CFT, mesons, baryons, finite density, confinement. 



* yseo@sogang . ac . kr 
Tshock@fpaxpl.usc.es 
■fsjsin@hanyang.ac.kr 
§ zoakos@fp axpl.usc.es 



Contents 



1. Introduction 1 

2. Baryons and baryon density in holographic QCD 3 

2.1 The confining geometry in holographic QCD 3 

2.2 Flavor in the confining geometry 5 

2.3 Quark density and chiral symmetry 7 

2.4 Baryon vertex 9 

2.5 Baryon density in a confining geometry 11 

3. Baryons in the dilaton flow geometry 12 

3.1 D5-brane setup 12 

3.2 Probe D7-brane 16 

4. Meson Spectrum 21 

4.1 Setup 21 

4.2 Zero m q spectral flows 24 

4.3 Large m q spectral flows 25 

5. Discussion 29 
A. Alternative form of the Schrodinger potential 30 



1. Introduction 

Quantum chromodynamics (QCD), known to be the microscopic explanation for the 
nuclear force, has been studied for more than 30 years. Nevertheless, reliable and 
fast methods for treating its dynamics in the low energy regime are still lacking. 
Although lattice QCD is rapidly developing, there are problems which this technique 
is not well suited to, in particular real-time calculations and calculations at finite 
density due to the infamous sign problem. We therefore need new ideas which work 
at least qualitatively in the presence of a chemical potential. Recently, techniques 
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derived from string theory have made a remarkable connection between strongly 
interacting gauge theories and gravity in asymptotically Anti De Sitter geometies [1] . 
Therefore it is natural to ask whether AdS/CFT can shed light on QCD dynamics, 
especially in the case of a dense medium. Such questions in the absence of finite 
baryon density have been partially answered in the approach known as holographic 
QCD [2, 3, 4]. Such avenues continue to provide new insights and new possibilities 
in the realm of understanding real QCD. The holographic encoding of the baryon 
chemical potential was studied in [5, 6, 7]. The equations of state were analyzed 
in [6] and the spectral functions were calculated in [8]. While the prescription for 
encoding finite temperature has been clear from the early days of AdS/CFT, it is 
still not clear in all contexts about the geometry describing finite baryon density. 
This is especially true in the case of QCD-like theories as will become apparent in 
this work. 

The introduction of fundamental matter by means of probe branes was pioneered 
in [9]. Such an addition allows for the study of meson phenomenology and the 
response of the system to the baryon density, encoded in a non-zero gauge field 
configuration on the flavor brane, is seen both in the flavor brane embedding along 
with the change in the spectral function. 

In the deconfined phase it has been shown that when quark density is introduced 
into the gauge theory (as encoded in a finite density of fundamental strings on the 
flavor brane worldvolume) the probe brane must end on the black hole horizon [7]. 
However, as pointed out in [10] when a compact D-brane corresponding to a baryon 
vertex is present (which is not possible in the deconfined geometry), a more natural 
situation is for the flavor brane to attach to the wrapped baryon vertex brane. In 
the same paper, the density dependence of the baryon mass was examined in detail. 

In this paper, we consider the meson spectrum in the presence of finite baryon 
density in a confining gauge theory following the approach developed in [10]. We use 
the dual confining background at zero temperature found in [11] with the addition of 
flavor D7-branes and a D5-brane baryon vertex joined by a force balancing condition. 
We examine the simplest mode on the D7-brane which describes the goldstone boson 
in the absence of a current quark mass m q . We find that if the current quark mass 
m q is sufficiently large, the meson mass reduces as we increase the baryon density 
while it exhibits a more complicated behaviour for small m q , which is similar to the 
behaviour of a baryon mass in the medium. In all cases we find that the baryon 
density causes a splitting in the spectrum and an intricate spectral flow which can 
be understood by studying the Schrodinger potential. 

The rest of the paper consists of the following sections. In section 2 the problem 
is set up by reviewing baryons and baryon density in holographic QCD. In section 3 
the interaction between the meson and the baryonic medium is discussed. In section 
4 the flow of the meson spectrum as a function of the baryon density is shown. In 
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section 5 we provide a discussion of the results provided here and add possible lines 
for future work. 

2. Baryons and baryon density in holographic QCD 

2.1 The confining geometry in holographic QCD 

We provide here a quick review of the Gubser dilaton flow geometry, [11], which we 
use as the holographic dual of our confining background. Similar non-supersymmetric 
geometries can be found in [13] and [14]. 

The Gubser dilaton flow geometry (GDFG) is a non-supersymmetric deformation 
of AdS§ x S 5 which corresponds on the field theory side to adding a vev for the 
supersymmetry breaking term (tr F 2 ) 1 . This operator is sourced on the supergravity 
side by the dilaton and because of the broken conformal invariance, the running of 
this operator is seen in the non-trivial profile for the dilaton. In the UV limit we 
recover the pure AdS geometry but in the IR the theory runs to a strongly coupled (in 
g s ) limit at which point the dilaton blows up, leaving a naked singularity. By probing 
this geometry with a fundamental string stretched from the boundary, corresponding 
to a Wilson loop it can be shown that the geometry presents an area law for the energy 
of the Wilson loop and therefore exhibits confinement. A mass gap is also found by 
calculating perturbations of the dilaton 2 . The position of the curvature singularity 
provides an IR scale in the field theory which we can think of as Aqqd which sets 
all further scales in the problem. The solution has a single dimensionful parameter 
which dictates the radius of the singularity. We can set this unique dimensionful 
parameter on the field theory side to 1 thereby scaling all dimensionful parameters 
in units of Aqcd- One clear problem in constructing a realistic theory from such 
a geometry is that there is no hierarchy between the strong coupling scale and the 
supersymmetry breaking scale. However, for the purposes of this investigation we 
will simply utilise the fact that this theory has confinement and a mass gap, just as 
in QCD. 

We clearly cannot trust the geometry in the region around the singularity but for 
the probes we will use to study hadronic physics the singular region exhibits a repul- 
sive potential and we can exclude all solutions where branes end on the singularity. 
Moreover, according to the criteria of [15] and [16], this singularity is good. The 
analysis in [15] concludes that the criterion is that the time component of the metric 

1 It should be noted that pure JV = 4 SYM contains no such 50(6) preserving operators but a 
deformation where an 50(6) invariant mass term for the scalars has been added may allow for such 
a term 

2 A further problem is that there is no separation of scales between i?-charged and non-charged 
glucball excitations 
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should not blow up at the singularity. Using the Einstein frame for this analysis, it 
is easy to prove that the GDFG geometry fulfills this criterion. 

The GDFG is a solution of the type IIB supergravity equations of motion and 
the solution can be written as: 



while the five-form remains unaltered from the pure AdS solution. R is the AdS 
radius and ro is the position of the singularity which we set to 1 in the following. 

The field theory dual of this geometry clearly has some strong similarities to QCD 
although there are some important differences. We would like simply to exploit the 
fact that we have a supergravity dual to a confining gauge theory, which is somewhat 
more realistic than the hard-wall models (see [3] and references therein) in order to 
further explore QCD-like behaviour. Some of the questions we might ask of such a 
model are: 

• Can we add flavor to such a model (both quenched and eventually unquenched)? 

• What is the finite temperature description of such a gauge theory? 

• What is the gravity dual of finite baryon density in a confining gauge theory 
and what is its phenomenology? 

Section 2.2 reviews the addition of quenched fundamental matter in such a gauge 
theory and it will be clear that some qualitatively QCD-like phenomenology comes 
out of such a picture. The unquenched calculation has not been studied but would 
certainly be an interesting direction for future work (see [17] for recent work on 
studying unquenched flavor in QCD-like models). 

In [18] and separately in [19] the finite temperature counterparts to this system 
were studied, concluding that at finite temperature the only possible solution to the 
supergravity equations with a horizon had to have a trivial dilation, and therefore 
no vev for (trF 2 ). This calculation however does not allow for the possibility of a 
solution with a naked singularity in the supergravity limit. 

The last question above is precisely what we would like to tackle in the current 
work. In a non-confining geometry one can add finite baryon density to a system 
simply with the addition of more quarks than antiquarks. We shall deal with this 
in section 2.3. However, the problem is more complicated in a confining geometry 
where we have to introduce a new object - the baryon vertex - in order to accomplish 
this. Clearly the addition of free quarks does not make sense in such a context. 
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where A(r) and the dilaton are given by: 
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2.2 Flavor in the confining geometry 



For discussions of flavor in the quenched approximation in AdS/CFT we refer the 
reader to [9, 20, 21, 22]. 

The addition of flavor has been studied extensively in the dilaton flow geometry 
in [23, 24, 25]. Here we remind the reader of the most important phenomenological 
features of a such a setup. We will be able to contrast these with the behaviour of 
the theory in the presence of a baryon vertex in the forthcoming sections. 

For convenience we write the metric (2.1) in a form where the natural embedding 
of a D7-brane is clearest. The radial direction and five-sphere combine to form an 
E 6 (up to confer mal scalings), written as E 4 x E 2 . The D7-brane lives perpendicular 
to the E 2 leading to a global U(1)r symmetry on the worldvolume in the case of 
massless flavors. The breaking of this U(l) by non-trivial embeddings of the D7- 
brane corresponds on the field theory side to the breaking of the chiral symmetry 
through the dynamical generation of a (qq) vacuum expectation value. 

We rewrite the metric as: 



ds -j^q — & ^ 



r 2 R 2 

—A 2 (r)r]^dx' 1 dx u + — (dp 2 + p 2 dti\ + dL 2 + L 2 dp> 2 ) 
R 2 r z x ' 



(2.3) 



where the E 4 and E 2 are manifest, r 2 = p 2 + L 2 and L and <p are the transverse 
directions to the D7-brane. Because of the manifest 50(2) symmetry of the geometry 
we are free to chose a solution with p> = 0. The induced metric on a D7-brane with 
Euclidean signature is: 



^A 2 {r){dt 2 + dx 2 ) + ^[{l + L 2 )dp 2 + p 2 dn 2 ^ , (2.4) 



where L = dL/dp. In order to study finite baryon density we will be interested in 
turning on a non- vanishing gauge field of the form Ao(p). 

The DBI action for the Nf D7-branes is: 

S D7 = -iV> 7 r 4 J dt s e~*> l /det( y g + 2Tia>F) 

= -NfPr J ^pA(r) 3 p 3 e* /2 y / e*A(r) 2 (l + L 2 ) - F 2 , (2.5) 

where we have denoted: 

r 7 = fi 7 Q 3 V 3 , F = 2na'F pt , (2.6) 

and we have turned the radial variables p and L into their dimensionless counterparts 
via (p, L) —>■ (r p, r Q L). We do not relabel these variables with a tilde as is often the 
case, for notational simplicity but in the future all dimensionful quantities originating 
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from the holographic energy/ radius duality will be in units of r when not explicitely 
stated. For the current illustration we will set the gauge field to zero. How to 
consistently turn this field on in a confining geometry will be the task of the bulk 
of this paper. At the classical embedding level we have a single equation of motion 
for the field L. In the large p, UV limit the solution has two free parameters, one 
normalisable and one non-normalisable: 

L(p) p ^oo = m + 4 • (2-7) 
P 

These correspond to the source and vacuum expectation value for (qq) (the quark 
mass m q = rom/27ra' and bilinear condensate (qq) ~ N c XtqC 3 ). Although there are 
two free parameters, the dynamics of the theory in the IR only allows well behaved 
solutions for a finite number of values of c for a given m, in the present case there is 
just one physical value of c for any m. In practice we work backwards and integrate 
out from the IR with the well behaved solutions given by L(0) = L c , L'(0) = 0. At 
the boundary we read off the values of m and c. These are parameterised by the one 
parameter family of solutions given by the value L c . 

In the following figures we show the embeddings for a range of L c values and 
the plot of the condensate versus the mass calculated parametrically in L c . We see 



L 

T- 




P 



Figure 1: D7-brane embeddings with different values of the IR boundary condition L(p = 
0) = L c . In this and all future such plots, the grey disk labels the singular region of the 
geometry. 

that there is spontaneous chiral symmetry breaking and therefore one may expect a 
goldstone mode. This can be found by studying the excitations of the field in the 
direction perpendicular to the condensate (as expected from goldstone's theorem). 
In the case of the D7-brane there are two scalar modes and an 8-compontent gauge 
field (before gauge fixing) which can be excited on the worldvolume. On top of 
the classical embedding we can study excitations which correspond to a higgsdike 
mode, L=L em b + 2ira'L. The goldstone mode is given by excitations around the 

3 for a description of the supersymmetrically complete operator corresponding to c see [7] 
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Figure 2: m versus c indicating spontaneous chiral symmetry breaking in the massless 
limit. 



classical solution which in this case has been set to zero by the underlying symmetry: 
if = + 2ira'<f). We will ignore fields with indices in the three-sphere directions, 
which are dual to R-charge currents, and set the gauge field compontent A p to 
zero as a gauge fixing condition. This leaves us with four vector fields. In the 
case that we study excitations with zero momentum in the Minkowski space, the S 3 
rotational symmetry leaves us with only two distinct gauge field components, A Q and 
Ai, i = 1..3. For the moment the interest lies in the Goldstone mode and we remind 
the reader of the result from [24] for the behaviour of the goldstone boson mass as 
a function of the quark mass. This follows the Gell-Mann-Oakes-Renner relation 
which can be shown exactly by using an expansion about the massless embedding 
(see [26, 27] for details). 

2.3 Quark density and chiral symmetry 

Now we turn to the study of finite quark/baryon density. This has been studied 
extensively both in the case of supersymmetric M = 4 SYM [28] plus flavor and the 
finite temperature counterpart [5, 6, 7]. 

Finite baryon density corresponds to a non-zero expectation value for the opera- 
tor g7°g which is sourced by the time component of the gauge field on the D7-brane. 
The solution to such a field Aq will generically have a UV behaviour containing both 
a normalisable and a non-normalisable piece. The non-normalisable piece is the 
source for the above operator which corresponds to the chemical potential, whilst 
the coefficient of the normalisable term is related to the baryon density itself. 

For illustrative purposes we will first discuss the addition of finite baryon density 
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in the case where there is a horizon in our space (including both the zero temperature, 
extremal case and the finite temperature non- extremal case). At this point it is 
perhaps more sensible to talk of a finite quark density which is related to the baryon 
density simply by a factor of the number of colors. On the D7-brane the non-trivial 
gauge field configuration acts to deform the brane, especially in the IR and the 
deformation manifests itself as a throat on the brane which ends on the horizon. By 
studying the tension of this throat it can be shown that it corresponds to a bundle of 
fundamental strings, pulling the D7- brane into the horizon, where the charges can be 
swallowed. Such a fundamental density of strings is expected due to the finite quark 
density, the charges of which have to end somewhere and not on the vanishing three 
sphere which the D7-brane otherwise wraps. Thus the brane behaviour is completely 
consistent with the field theory picture. 

The question of interest in the current paper is what happens to a confining 
gauge theory when finite quark density is added. In the dilaton flow geometry dis- 
cussed in the previous sections, a D7-brane embedding which falls into the naked 
singularity would clearly not be a good solution. We can then ask if there is any 
other configuration which is well defined in the presence of finite quark density. 

For the case of pure AdS 5 x S 5 we illustrate the embeddings studied in [28]. In 
this case the action for a D7-brane with finite baryon density is given by: 

S D7 = -t 7 J dtdpp 3 ^(l + L 2 + LV) - A> ( P y . (2.8) 

Again we use the U(l) symmetry to set ip — 0. We also note that there are two 
conserved charges as both L and Aq appear only with derivatives. The gauge field 
A corresponds, as described above, to a finite baryon density. The conserved charge 
with respect to the gauge field is given by: 

Q=- . = , (2.9) 



Performing the Legendre transform and replacing the gauge field with the conserved 
charge we are left with: 



S 



D7 



-r 7 J dtrf PV /(l + L 2 )(Q 2 + p 6 ). (2.10) 



Note that there is a trivial scaling where we can rescale all dimensionful quantities 
either in units of Q or in units of the asymptotic value of L which gives the quark mass 
(or the non-normalisable behaviour giving the condensate). The physically relevant 
parameter we can tune is then Solving the equations of motion numerically and 
plotting the D7-brane embeddings in figure 3 we see a family of solutions which have 
the clear throat-like behaviour in the IR. In the case of finite temperature there is 
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an extra scale and so the possible phase space is more interesting, though the same 
general tendency holds that the fundamental strings dissolved on the D7-brane pull 
it into the horizon, giving the throat-like behaviour. 



L 




p 



Figure 3: D7-brane embeddings in the case of pure AdS§ x S 5 with finite baryon density. 
In this case the only dimensionless scale in the problem is the ratio of the baryon density 
to the quark mass. 



2.4 Baryon vertex 

Having introduced finite baryon density in the previous section and realised that 
there may be a problem in the case of a confining geometry, we now turn to the 
study of the baryon vertex, first introduced by Witten in [29]. The baryon vertex 
is a gauge invariant antisymmetric combination of N external quarks, whose dual 
gravity description is given by a D5-brane wrapping a five-sphere at some radius in 
the AdS space and connecting with the boundary by N fundamental strings. (NB. 
we stick here to the case of AdS/CFT in the AdS§ x S 5 context) 

Witten's argument goes as follows: In type IIB string theory there is a self-dual 
field strength F 5 and the compactification on AdS$ x S 5 provides N c units of flux on 
the five-sphere, J g5 = N c . On the D5-brane world volume there is a U(l) gauge 
field A which couples to the five-form field strength through the term J RxS5 A A It 
is because of this coupling that G 5 contributes N c units of U(l) charge. Each string 
endpoint adds —1 unit of charge and since in a compact space the total charge has 
to vanish, N c strings have to end on the D5 brane. In the SU(N C ) gauge theory the 
gauge invariant combination of N c quarks is completely antisymmetric and, indeed, 
the strings between the boundary (or a D3-brane) and the D5-brane are fermionic 
strings. 
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There are two complementary approaches to the baryon vertex depending on the 
influence of the N c fundamental strings on the D5-brane. In the first one we neglect 
the deformation of the D5-brane and the world- volume gauge field on it, due to a 
uniform distribution of fundamental strings over the five sphere. In this way we treat 
the configuration as a combination of strings with a flat D-brane wrapping around 
the S 5 , for the purpose of calculating the total energy of the brane [30, 31]. However, 
the distribution of N c fundamental strings breaks supersymmetry completely. Even if 
each string preserves one half of the supersymmetry, the fact that they take different 
positions on the D5-brane makes the preserved Killing spinor of each string different, 
such that all supersymmetry is broken [32]. 

The supersymmetric solution is to have the D5-brane deformed by means of the 
tension of the strings attached to it. This deformation will happen if a significant 
number of them join at the same point. In this way the string back-reaction is 
not negligible and should be taken into account [31, 32]. In this case the full DBI 
action should be considered, where the fundamental strings are seen as a spike in the 
world- volume of the D5-brane, along the lines of [33, 34]. 

In the approaches already discussed the quarks are non-dynamical and are de- 
fined by strings stretching from the D5-brane to the AdS boundary. In this way they 
are describing infinitely massive quarks. The question then is whether one can add 
dynamical (though not in the sense of unquenched) quarks with the usual probe em- 
beddings of D7-branes and have the fundamental strings from the D5-brane attach 
to the D7's. On the gauge theory side (in the pure AdS limit) we will have M = 2 
supersymmetry and fundamental quarks, so we expect dynamical, finite-energy, su- 
persymmetric baryons. The gravity dual of such an object will be a D5-brane wrap- 
ping a five-sphere and connecting to the D7-brane with iV c fundamental strings, [20] . 
The fact that the baryon vertex is supersymmetric (in the pure AdS background) is 
related to the orientation of the D7-brane. As long as it is parallel to the D3-branes 
and orthogonal to the strings the dynamical baryon will preserve half of the ini- 
tial supersymmetry. The situation changes dramatically upon introducing a nonzero 
temperature. In this case, the gravitational attraction of the black hole pulls the 
D5-brane towards the event horizon until it collapses, leaving us with just a black 
hole embedding for the D7-branes, [38]. This is to be expected as we already know 
that when we have a finite density of fundamental quarks dissolved in the brane it is 
pulled into the black hole under the tension of the strings, so the situation with and 
without the baryon vertex at finite temperature is identical. This is not surprising 
from the field theory point of view as in the deconfined phase we don't expect to 
have any baryons. 

In the Sakai Sugimoto model the picture is very similar to the situation described 
above with a wrapped D4-brane mirroring the wrapped D5-brane in our description 
[37]. In this context the holographic dual to the skyrme model has been constructed 
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whereby the baryon is built as a topological configuration of the pion degrees of 
freedom. The spectrum of baryons is then found by solving the Schrodinger equation 
in the potential which defines the moduli space of the skyrmion. An analogous picture 
has not been constructed in the current setup. 

2.5 Baryon density in a confining geometry 

At this stage we can motivate our setup before going onto the full construction in 
the following sections. 

We know that if we add a finite baryon density to a D7-probe brane, this cor- 
responds to having a density of fundamental strings dissolved in the brane world 
volume. The charges from these have to end somewhere and in the case of finite 
temperature the natural place for the strings to end is on the horizon. The baryon 
vertex can be thought of as a density of fundamental strings dissolved on the world- 
volume of a wrapped D5-brane and these also need to end somewhere. In the trivial 
setup these end at the AdS boundary corresponding to a bunch of infinitely massive 
fundamental strings. However, if we want to introduce true quarks, corresponding to 
dynamical objects in the fundamental of a global symmetry and the fundamental of 
the color group these strings should end on a flavor brane. In a confining geometry 
it seems natural to connect the fundamental strings from the D5-brane with those 
from the D7-brane. 

We can now motivate the computations which follow by a simple cartoon as 
shown in figure 4. The image on the left illustrates the most trivial, non-realistic 
setup. A D7-brane brane with finite charge density looks to have a throat made 
of fundamental strings which wants to end somewhere to carry off the charge. A 
baryon vertex made of a D5-brane wrapped on the five-sphere also has a throat 
made of a bunch of fundamental strings which wants to end on the AdS boundary. 
The natural interpolation (given that the fundamental strings pull the two objects 
together) is then the right diagram which shows the D5 and D7-brane meeting at 
a point, with an appropriate force balancing condition to create a stable situation. 
However, this may not be the true lowest energy solution but simply a metastable 
configuration. It seems likely that given the point-like nature of the vertex there 
may be some non-trivial correction to this which will be the true vacuum. The case 
at hand is different from that of the Sakai-Sugimoto DA — D8 — D8 case as the 
D4-brane vertex is completely within the D8-brane worldvolume and therefore an 
instanton interpretation is natural. It would be interesting to know if there is an 
analogue to the skyrme model discussed in [37] in the present context. See section 6 
in [20] for a discussion of the difficulties involved in such a calculation for the present 
case. 
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Figure 4: Left image: Blue D7-brane with a bunch of fundamental strings leading off it in 
the form of a throat plus red, wrapped D5-brane with a similar throat. Right image: The 
two throats meet at a point in the (p, L) plane (the only allowed configuration in terms of 
the codimension of the objects). The configuration should be force balanced between the 
two objects. 



3. Baryons in the dilaton flow geometry 
3.1 D5-brane setup 

Here we study the embedding of a wrapped D5-brane in the GDFG. This will be 
the baryon vertex which, under force balancing conditions will attach at a point 
to the D7-brane, providing us with a baryon vertex where the constituents can be 
understood as dynamical quarks. 

The background is given by the metric (2.1) and we will set to 1 in all further 
numerical calculations. For the wrapped configuration of the D5-brane, the back- 
ground five-form field strength (sourced by the D3-branes) will couple to the world 
volume gauge field Aq via a Wess-Zumino term on the D5-brane worldvolume. The 
five-sphere coordinates are split into 9 and Q4 and we will find solutions where the 
scalar field on the D5-brane corresponding to the radial direction of AdS is a function 
of 9 only. The gauge field will also be a function of this direction only. This leaves us 
with an SO(5) symmetry for the wrapped, deformed D5-brane. We therefore have 
r = r(9), Aq = Aq(9). Such a setup corresponds to having the bunch of fundamental 
strings attached at a single point on the D5-brane world volume. The r direction 
corresponds in the D7-brane language to ^/p 2 + L 2 . In the end the two solutions 
will join up at a point where p = meaning that L = r at the vertex. 
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The induced, string frame, metric on the D5-brane with Euclidean signature is: 
~ 2 R 2 



d? 2 

ab D5 



A{r) 2 dt 2 



R 2 



( r ' 2 + r 2 ) d9 2 + R 2 sin 2 9d£l 



(3.1) 



where r' = dr/d6. The DBI action for single D5-brane with N c fundamental strings 
can be written as: 

S D5 = -n 5 J d^e~' s ' v / det(g + 2TTa'F) + fi 5 J 2na'A A G (5) 

i : : 

(3.2) 



7 -, / IdO sin 4 9 -\Je*A(r) 2 (r 12 + r 2 ) — F 2 + AA 
= r 5 J dtCob , 

where we have denoted: 

r 5 = /i5 fi 4j R 4 , A = 2naA k F — A' {r) . 
The dimensionless displacement can be defined as follows: 

dC D5 sin 4 OF 



(3.3) 



dF 



e*A(r) 2 (r' 2 + r 2 ) - F 2 



= -D(9) , 

and the equation of motion for gauge field is given by: 

d e D(6) = -4sin 4 fl. 

Integrating, we have: 

D{9) 



3 3 

- [y-K — 9) + - sin 9 cos 9 + sin 3 9 cos 9 , 



(3.4) 



(3.5) 



(3.6) 



where the integration constant v determines the number of fundamental strings [y N c 
strings are attached to the south pole and (1 — v) N c strings to north pole). This 
solution holds true in any geometry where the E 6 in the metric is not deformed (up 
to conformal rescalings). By performing a Legendre transformation with respect to 
the gauge field, we arrive at the Hamiltonian: 



l~lD5/ r 5 — r *-> D5 

OF 



= J d9A{r)^e^{r' 2 + r 2 )^D{9) 2 + sin 8 9, 
the equation of motion for which is: 



(3.7) 



d 
d9 



r'A(r)e^ 2 y/D(9) 2 + sin 8 9 
\Jr' 2 + r 2 



d_ 
dr 



A(r)e* /2 v / r /2 + r 2 yj D{0) 2 + sin 8 9^j =0. 

(3.8) 
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There is a trivial solution to this equation of motion given by r = ~ 1.471 in 
units of r . This is the first artifact of the confining geometry. In pure AdS or 
AdS/Schwarschild there is no stable solution for a wrapped D5-brane as the mini- 
mum energy configuration is always when the brane collapses into the horizon (be it 
extremal or otherwise). 

In the following we set v = 0, meaning that all fundamental strings attach to 
one pole of the D5-brane, 9 = ir. The equations of motion can be solved very simply 
using any numerical integration package. The boundary conditions which we impose 
are set at 9 = 0. These are labelled = r{9 = 0), setting also r'{9 = 0) = 0. Solving 
with these boundary conditions we find a range of different behaviours depending 
on the value of ?y For small values of the D5-brane falls into the singularity, and 
clearly such solutions should not be trusted in the supergravity limit. For larger 
values of but less than r* the solutions wrap around the singularity and end at 
9 = 7r at a value labelled r c with a positive gradient, r'{9 = n) < 0, whereas for 
solutions with r« > r* the solutions end at r c with r'(9 = ir) > 0. This change in 
behaviour will be important in future sections. In figure 5 we plot a range of different 
solutions including the singular solutions, and the 'well-behaved' solutions including, 
in black, the r* solution. 




Figure 5: Wrapped D5-brane solutions for different rj. The black circle corresponds to 
the r* solution. 
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The hamiltonian (3.7) gives the energy of the D5-brane. The r c dependence of 
the mass of a single D5-brane is drawn in figure 6. We see that as the value of r c 
goes to one, the mass increases but does not go to infinity as the brane touches the 
singularity (though the Hamiltonian density does diverge on the singular surface but 
slower than the worldvolume of the D5 vanishes). It should also be noted that in the 
absence of any external forces, the equilibrium position for the brane is at r = r*. 

M D5 




Figure 6: r c dependence of the energy of the wrapped D5-brane. The minimum gives the 
value r* for which the pull of the singularity exactly cancels the pull of the fundamental 
degrees of freedom at the North pole of the brane. 



Because any D5-brane sitting away from r = r* is out of equilibrium there is 
a force at the apex r c where the tension from the brane tries to pull it back to its 
equilibrium position. This force can be obtained from the variation of Hamiltonian 
of D5-brane with respect to r c : 



'D5 



dr c 

A(r)r'e $ / 2 , , 

(3.9) 



-NrT, 



r=r c 



where Tp is the tension of a fundamental string and the direction is towards the 
singularity at the centre of the geometry (though note that for solutions with negative 
gradient r c < r* the force will be in the opposite direction). In order to match this 
force we can place another object which also has N c fundamental strings attached 
to it pulling in the opposite sense [10]. This is quite a natural set-up. The D5-brane 
has a bunch of fundamental strings attached in order to carry away the flux flowing 
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through its world volume. We can let these strings go to the boundary in which case 
we can consider the system a baryon with infinitely massive quarks, or we can let 
them end on a flavor brane. The flavor brane which in this case is a D7-brane can 
have fundamental strings ending on it if we turn on a finite baryon density as we 
have discussed above. Eventually the stable equilibrium of the system will be one 
where the fundamental strings have shrunk to zero size, pulling both the D5 and the 
D7-brane to a cusp. 

3.2 Probe D7-brane 

Now, we reconsider the D7-brane probe brane with finite baryon density in the 
confining geometry with the appropriate force balancing condition. The D7-brane 
solution will have a cusp where the fundamental strings from the D5-brane attach. 
These end points can be understood as point charges on the D7-brane. As before, 
the action for Nf D7-branes is: 

= -t 7 J dtdpA(r) 3 p 3 e^ 2 ^e"A(r) 2 (l + L 2 ) - F 2 

= t 7 f dt£ m , (3.10) 



where we have denoted: 

r 7 = p 7 Q 3 V 3 & F = 2na'F p0 . (3.11) 

Just as in the pure AdS geometry we can look at the equation of motion for the 
gauge field to calculate the dimensionless quantity Q: 

8C m _ A(r)Ve^P = -QT F = -Q, (3.12) 



d ^ \J e®A{r) 2 {\ + L 2 ) — F 2 
which can be replaced in the Legendre transformed hamiltonian to give us: 

= J dp^Q 2 + A(r) 6 p 6 e*y/e*A{r) 2 (l + L 2 ) . (3.13) 

Just as before we can solve the equations of motion numerically, however now we 
have an extra force on the D7-brane coming from the fundamental strings stretching 
between it and the D5-brane. The force at the cusp of the D7-brane is equal to: 

dH D7 
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_ QT Mr)Le*l> 



(3.14) 



L=L C 



where L c is the location of the cusp. This must be balanced with the force from the 
D5 and so we get the condition that: 

F D7 = ^F D5 , (3.15) 

This states that the force on the D7-brane increases as we increase the quark density 
{n q ~ ^2-). The above relation can be written in the following simple form: 



(3.16) 



where we have used the fact that r c = L c . Now, solving the equations of motion with 
this as a boundary condition in the IR we can calculate the D7-brane embeddings 
for a given value of the quark mass and baryon density. In practice the calculation 
works in a different order: we pick a value of and solve the D5-brane embedding 
to find r c (rj). Using this value, along with the behaviour of the D5 near the cusp as 
the boundary conditions for the D7-brane with a given baryon density we integrate 
this out to the boundary and read off the value of the quark mass. We implement 
an efficient sequence of binary search routines to find an embedding with a fixed 
value of the quark mass as we vary Q. Now we can get various D7-brane embedding 
solutions with different choices of parameters. 

In figure 7 we plot a range of m ? 's for different values of Q. At each cusp at 
p = there is a balancing of forces between the D5-brane and D7-brane. It can be 
seen by comparing the two graphs in figure 7 that larger Q affects the behaviour of 
the D7-brane from a larger value of p whereas for small Q the effect is only seen at 
small p. It is interesting to see how we must change r^, the South pole position of 
the D5-brane as a function of Q for a fixed mass. In figure 8 we plot this for the 
massless embeddings and note that there is a nearly logarithmic dependence on Q. 
Such a dependence will be seen in many plots which follow. 

In figure 9 we plot a series of m q = embeddings with Q ranging from e -6 to e 6 . 
We see that some of the D7-brane embeddings intersect with the D5-brane at a point 
other than the apex. Such solutions should perhaps not be trusted as there ought to 
be some non-trivial interaction between the two branes at these intersection points. 
The question of whether or not the D5-brane intersects the D7-brane is a simple one 
to answer. Because of the force balancing condition the gradient of the D7-brane 
at the cusp is related up to a scaling, to the gradient of the D5-brane at the cusp. 
Only when the D5-brane has a positive curvature (in the coordinate system drawn 
in the figures here) will there be no intersection. The point at which this change 
occurs is at r = r* and therefore we can parametrically plot the graph (m g (r*, Q), Q) 
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Figure 7: m q = 0, ...,3 with (a) D7 and D5-brane embeddings for Q=0.1 (b) D7 and 
D5-brane embeddings for Q=5. 



ri 




Log[Q] 



Figure 8: Q dependence of the position r, of the South pole of the D5-brane. We see 
there is very slow dependence of r, on Q, however, as can be seen from figure 4, r c depends 
very strongly on the position of r^. 

(figure 10) which defines the line separating the region of intersecting solutions from 
the non-intersecting solutions. It seems likely that the intersecting solutions are less 
physically realistic and so we will concentrate on the non-intersecting regions where 
possible. We do note however that the Q — > limit recovers smoothly all known 
solutions in the absence of the baryon vertex and baryon density. 



- 18 - 



L 

2 r 




Figure 9: m q = plots for a range of LogfQ], from lowest D7-brane to highest Log[Q] = 
(—oo,—5.4,—3.4,—1.4,0.6,2.6,4.6,6.8) - D7-branes are blue, D5-branes are red. We see 
that only for large values of Q in this case > 10 do the D5 and D7-branes not intersect 
each other away from the apex. 




Figure 10: Region of intersections of the D5-brane and D7-brane shown in the blue shaded 
area below the line. The line corresponds to the embeddings with rj = r*. 
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In figure 11 we plot another set of D7-brane embeddings, but this time for fixed 
values of m q in order to see the change in behaviour with Q. It is shown that for 
small and large quark masses, the effect of adding baryon density is different. For 
small m q turning on Q increases the 'dynamical quark mass' (increases the value L 
at which the D7 and D5-branes meet), whilst for large quark mass the 'dynamical 
quark mass' is decreased. Note that the turnaround point where increasing Q starts 
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Figure 11: top left to bottom right: D5/D7-brane solutions for m q = (0,1,2,5) for Q 
in the range Log[Q] = (— oo, —5.4, —3.4, —1.4, 0.6, 2.6, 4.6, 6.8). The D7-brane embeddings 
are color coded starting with blue for small Q and going to pink for large Q. It can be seen 
that for smaller m q , increasing baryon density lifts the value of L c whilst the opposite is 
true for larger m q . 



to drag the D7-brane down rather than up in the IR (for a given m q ) is around 
m q = 1.6. 

For each of the sets of embeddings we can study the asymptotic behaviour of 
the D7-brane solution and, just as in the Q = case we can read off the value of the 
condensate. In figure 12 we plot the condensate as a function of the mass for two 
different values of Q. We note that there is a chiral condensate for zero quark mass. 



Figure 13 shows the Q dependence of the energy of the D5-brane from the integral 
of the hamiltonian density for fixed m q . We see in figure 13 that for very small values 
of the quark mass there is a minimum value of the mass of the D5-brane which is 
absent in the D4/D6 system [10]. 
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Figure 12: Chiral condensation for Q = 0.1 (blue line) and for Q = 5 (red line). Given 
that at finite baryon density in a non-confining geometry a chiral condensate with the 
opposite sign to that in a confining geometry at zero baryon density appears, a cross-over 
in behaviour is not surprising. 



Figure 13: Density dependence of the D5-brane energy for different m q . The left figure 
is (from bottom to top at log(Q) = 7) m q = (0.1, 0.3, 0.5, 0.7, 0.9, 1.2, 1.6) and on the right 
m q = (1.8, 2.5, 3.5, 4.5, 6, 8, 10, 13, 17). Again the very slow dependence on Q is clear here. 

4. Meson Spectrum 
4.1 Setup 

At this point we turn our interest to studying the spectrum of mesons in our setup in 
the presence of a baryon vertex and baryon density. Having calculated the classical 
solutions we can ask about the excitations on top of the embeddings. The Minkowski- 
spacetime dependent excitations which have normalisable boundary behaviour in 
the UV correspond to excitations of the quark bilinear operators on top of any 
possible condensate value. These therefore correspond to mesons. Depending on the 
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field we chose to excite on the D7-brane we will be able to study scalar and vector 
representations of mesons. The excitations are treated as small perturbations (of 
order a' ~ and therefore we may expand the action up to quadratic order in these 
modes. Upon imposing the equations of motion for the background fields, the action 
linear in fluctuations vanishes and the equations of motion for the perturbations can 
be derived from the quadratic piece only. 

The fields which may be of interest in the subsequent analysis are the gauge field 
and the two scalars (corresponding to excitations of scalar and pseudoscalar quark 
bilinears). We make the gauge choice of A p = and write our solution ansatz such 
that they depend only on p and x$. We label the fluctuations Ad,i,2,3(p, xq), L(p, 2 ) 
and (p(p,x ). Clearly one could study the system where the fluctuations have finite 
momentum within the finite density medium, but this will unnecessarilly complicate 
the current analysis and so for the time being we concentrate only on finite frequency, 
zero momentum solutions. 

The finite baryon density causes a coupling between the time component of the 
gauge field and the L fluctuations. We are however most interested in the pseu- 
doscalar which corresponds to the goldstone mode in the Q = case and this 
mode decouples completely even at finite baryon density. The action for this mode 
at quadratic order is given by: 

A(0) = \ [f( P ) d p 0(p, x ) 2 - h(p) d X0 0(p, x ) 2 ] , (4.1) 
where we have denoted: 



Hp) - A{pW + L {p?? V (Q 2 + p 6 Mp)^)(i + L{ P Y) . (4.2) 

Calculating the equations of motion we use the following parametrization for the 
scalar: 

4>(p, x ) = (p(p) e~ iwxa with u? = M 2 , (4.3) 
and we arrive at the equation for the (p(p): 

d p [f(p)dMp)} + M 2 h(pMp) = . (4.4) 

It will be convenient to calculate the solution of the equation of motion after we 
transform it to a Schrodinger form, where we have found that the analysis is nu- 
merically more stable than in the canonical form 4 . In order to convert (4.4) to a 



4 We will in fact use two Schrodinger forms in what follows, one of which is simpler for performing 
the numerical calculations and one of which gives a clearer picture of the bounded potential. The 
latter is introduced in the appendix. 
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Schrodinger form we decompose (p(p) as: 

<p(p) = ' (4 - 5) 

and arrive at the following differential equation for ip{p): 

iP" -V ef fip = 0, (4.6) 

where the effective potential is given by: 

1 f" 1 f 2 h 

Expression (4.7), contains terms which have second and third derivatives in the 
background embedding. Such derivative terms are numerically very unreliable, so 
we remove them using the equation of motion for the embedding to write it purely 
in terms of zeroth and first derivatives. 

In order to study the meson spectrum we study the excitations of the normalis- 
able modes. Doing so we perform an asymptotic expansion of the equations in the 
UV and find that the potential and the corresponding solution in this limit behave 
as: 

V p tL~lj 2 & ^-oo~Apt + -^. (4.8) 

We therefore take the normalisable solution and set the UV boundary condition 
at some puv to be ip(puv) — {Puv)~ 1 ^ 2 - I n the IR we are supposed also to make 
sure that our wavefunction is normalisable. Since the potential goes to a negative 
constant, depending on the position of the D7-brane and its derivative, the solution 
behaves as: 

V p e l{ ~ C & V^o ~ a cos WCp\ + c 2 sin U/Cp] , (4.9) 

c\ and C2 being integration constants. Calculating the integral of the square of 
the wave-function over the volume form it can be shown that both solutions give 
convergent IR behaviour and so are normalisable, meaning that we will have separate 
wave-functions with both Neumann (TV) and Dirichlet (D) boundary conditions. In 
the case of D7-branes in pure AdS, N and D boundary conditions coincide. The IR 
behavior in pure AdS is given by a sum of positive and negative exponential factors, 
so the normalisable solution has both vanishing value and derivative. In the current 
non-AdS case with finite Q, because of the two distinct boundary behaviours we may 
have a splitting in the spectrum (different spectra for N and D boundary conditions). 
This will correspond to a splitting in parity modes for the + and 0~ mesons. In 
the case of zero baryon density where the wavefunction for the meson mode is even 
about p = the mode is pseudoscalar because of its transformation properties in the 
(p, 0) plane (see [26] for a detailed discussion). 
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4.2 Zero m q spectral flows 

We now study how the baryon density in the confining geometry affects the goldstone 
spectrum. In order to do this we study the m q = solutions as a function of Q. In 
the Q = case we have a massless pseudoscalar mode which follows a Gell-Mann- 
Oakes-Renner relation in the small m q limit. When we turn on a finite baryon density 
we find the spectrum shown in figure 14. There are several key features to this plot. 




Log[Q] 



Figure 14: m q = flows as a function of Q. We see clearly the splitting between the 
Neumann and Dirichlet modes and the tachyon in the lowest Neumann mode as the eigen- 
value crosses zero (into the complex M n plane). (Neumann in red, Dirichlet in blue). The 
Neumann flow in the right figure continues to smaller Q (tending to M n — > 0) but is not 
shown in this graph. The right plot is a zoom in on the small Q region of the right plot. 

The first bulk feature is that there is a splitting in the spectrum between N and D 
modes which coincide at Q = 0. For large Q the modes split and there seems to be 
an equal splitting between each mode. 

There is also a clear difference in behaviour between iV and D modes for small 
Q. The most important difference is that the N modes contain a tachyon, indicating 
a clear instability above some value around Q = 2. This tachyon, being in the 
q^q part of the spectrum seems to indicate that the system wants to condense this 
operator. This would correspond to breaking of parity invariance (see [35] for details 
about parity conservation in vectordike theories). However, it is believed that in 
QCD there might be such a parity breaking phase transition at high baryon density, 
it would be extremely interesting if we are seeing the signature of such a breaking 
in this context [36]. It seems that we are getting a splitting in parity pairs by the 
presence of the baryon vertex. Indeed the baryon vertex gives an explicit breaking 
of the chiral symmetry (therefore lifting the Goldstone mode) and allows for the 
splitting of the + and 0~ modes. 

In fact the value of Q for which we get the tachyon is about the same value of 
Q where the intersections between the D5 and D7-brane disappear as the gradient 
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of the D7-brane in the IR becomes positive definite. We can ask how the tachyon 
behaviour depends on the value of m q and this is plotted in figure 15. 




Figure 15: A plot showing the region with the tachyonic instability and the region with 
intersections of the D5 and D7-brane. The wedge in the middle is the only region which 
we believe to be stable from the current analysis. The regions with the intersection and 
the tachyon appear to coincide for m q = at Q ~ 2.). 



4.3 Large m q spectral flows 

We now turn to the limit of large m q . In this limit at zero baryon density the D7- 
brane is unaffected by the deformation from pure AdS and so the spectrum is known 
analytically to be M n = 2m q ^(n + l)(n + 2). We can now follow the deformation 
from this result as we turn on the baryon density. We choose a quark mass m q = 20 
which at Q = gives very good agreement with the pure AdS result. 

In figure 16 we plot this spectrum and notice both some differences and similar- 
ities from the m q = case. 

The first similarity with the m q = plot is that there is a breaking in the 
degeneracy between N and D modes. The clear difference is that the general trend 
of the spectrum is downwards in contrast to the m q = case. However, this can 
be resolved by seeing the trend in the dynamical mass of the quarks as a function 
of Q for large and small m q in figure 11. It should also be noted that there is a 
tachyon appearing at around e 9 . This number seems very large, but because the 
units of Q are energy 3 and the scale of m q is 20 in the current setup, the value of 
e 9 ~ 8000 ~ 20 3 is not unreasonable (Note that for small m q the other important 
scale in the problem is ro and so figure 15 doesn't exhibit a cubic scaling for small 
m q . Fitting the large Q behaviour we find that roughly M n ~ ^ 2 ^.) 
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Figure 16: Spectral flow of the first six eigenstates as a function of Log(Q). We see that 
the AdS solutions (given by dotted lines) act as attractors for the spectral flow, all the way 
to very large baryon density. As the value of a meson mass gets close to that of the pure 
AdS case, the wavefunction is only slightly deformed about the AdS dip in the potential. 
As the flow goes between these values the wavefunction is affected more by the second 
potential well and finally for very large values of Q this is the dominant feature. Indeed 
we can continue this figure in the large Q direction and find that the first mode becomes 
tachyonic at around log Q = 9 - again it should be noted that this very slow dependence on 
Q is related to the dependence of the D5 embedding on rj. Again, red is Neumann, Blue 
is Dirichlet.). 



The other clear feature of the graph is the spectral flow whereby the spectrum 
moves between the pure AdS behaviour in steps. We can see this by looking at the 
Schrodinger potential in the form explained in the appendix in section A. Note that 
this flow is extremely reminiscent of the spectral flow on the Higgs branch studied 
in [39, 40]. 

In order to understand the behaviour we first study the Schrodinger potential for 
the case of pure AdS with zero baryon density. In this case there is a single minimum 
and the potential is given by the expression: 

Ae 2 ^m 2 (n + l)(n + 2) 

V n = l q - — 2 (4.10) 

{m 2 q + e 2 ^) 

This is plotted in figure 17. 

On introducing a non-zero baryon density, the D7-brane is bound to connect 
with the baryon vertex and so we get a radical change in the behaviour of the 
potential. The change in the embeddings can be seen in figure 18 illustrating the 
change in embeddings for a range of Q. We see that for small Q, the embeddings 
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Figure 17: AdS potential for the first six meson masses for m g = 20. 



are less deformed from the Q = embedding and the deformation only appears 
at small p. The effects of this will only be observable for high energy solutions 
whose wavefunctions have non-trivial support mostly in this region. Such high energy 
solutions are the higher resonances of M n . We therefore expect that as we increase 
Q from zero, the higher resonances will be affected first, just as we see in figures 16. 



Figure 18: D5-D7-brane bound solution for m q = 20 and log(Q) = (—6, —4, —2, 0, 2, 4, 6). 



The potential in the small Q limit goes like: 
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The IR limit of this has now changed from Vm = 1 to Vm = 1/4 and in accordance we 
have to alter the IR boundary conditions for any excitations. However, the dominant 
region is still the pure AdS potential well and so the spectrum remains the same as 
the AdS potential in the very small Q limit (except for high energy resonances). The 
new feature of a maximum in the deep IR only affects these high energy solutions. 
We plot the potential for very small baryon density in figure 19. As we increase 



v„ 




Figure 19: Schrodinger potential for the first six meson masses for an infinitesimally 
small baryon density. We see that the new feature is a very small perturbation on the AdS 
potential at around £ = —2. 

Q we see that the new maximum moves further into the UV, as the D7-brane is 
pulled more strongly towards the baryon vertex. We can see this from the brane 
embeddings, as shown in figure 18. 

The evolution of the potential with Q involves two factors. The first is that the 
new feature in the potential moves further into the UV, and the second is that a 
new minimum appears next to the maximum which slowly becomes an important 
feature in the potential. In figure 20 we plot the potentials for the first six modes for 
varying Q. If we plot the spectra along this flow in Q as we traverse the behaviour 
between the single potential and the double potential form we see that the meson 
masses are most concentrated around the AdS behaviour, stepping quickly between 
the AdS eigenstates. As we increase Q, although the meson masses are going down, 
there is a step like behaviour and the n th eigenstate takes the value of the n — m th 
AdS eigenstate in m steps until we lose completely the AdS-like behaviour for very 
high values of Q. The reason for this step like behaviour is because as the dynamical 
quark mass decreases (the D7-brane is pulled further into the IR) with increasing 
baryon density, the spectrum decreases. Every time it gets close to the AdS spectrum 



- 28 - 



Figure 20: variation of the potential as a function of Q. Each plot contains the potentials 
of the first six eigenvalues and the four plots run from top left to bottom right with 
log(Q) = (—4,0,4,8). We see how the two competing potential wells dominate in different 
ranges of Q values. (Neumann in red, Dirichlet in blue). 

it sits comfortably within the original AdS minimum and is not highly effected by 
the extra feature (the new minimum). 

5. Discussion 

In this paper, we have studied the meson spectrum in the presence of a finite baryon 
density in a confining gauge theory. The gravity dual of this gauge theory is a ge- 
ometry with a flowing dilaton and a naked singularity giving a scale for a mass gap. 
The interaction between the fundamental degrees of freedom and the baryon density 
is encoded in the interaction of the flavor brane and baryon vertex, modelled by a 
wrapped D5-brane. We have examined how the meson spectrum flows as we increase 
the baryon density and found that for sufficiently large current quark mass, the mass 
spectrum decreases as the baryon density goes up. This is seen on the flavor brane 
by a decrease in the dynamical quark mass. For near zero current quark mass, the 
spectral curve does not show monotonic behavior. One of the most interesting phe- 
nomena discovered in the current work is the presence of a critical baryon density 
where the meson mass vanishes and subsequently becomes unstable. The presence 
of such a density is robust regardless of the current quark mass, though its magni- 
tude depends on the value of m q . It would be interesting to study this further in 
the context of chiral symmetry restoration. Such an instability seems to signal the 
condensation of a pseudoscalar biquark operator though the true groundstate in this 
regime appears to be a more complicated interaction between the D7-brane and the 
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D5-brane than we have been able to model. A natural extension to the current work 
would be to look for the stable solution for the large baryon density regime. 

One of the new features in the calculation of the meson spectrum is that the IR 
boundary conditions for the D7-brane fluctuations are split into even and odd parity 
modes in the radial AdS direction. This is brought about by the breaking of the chiral 
symmetry explicitly by the baryon vertex. It is clearly difficult to phenomenologically 
match these results with those of real QCD as experimental results at finite baryon 
density but zero temperature are not available. Such a situation would occur in cold, 
dense objects such as neutron stars and so any spectral signatures from such objects 
would be extremely interesting. 

In addition to the phenomenology studied here one could think of using the cur- 
rent setup as a model for more, realistic QCD-like phenomena. Clearly in a realistic 
dual of QCD (beyond the large N c limit) we expect to see mesons with a finite width. 
One possible way to simulate such a width without looking at large N c corrections 
to the geometry, would be to allow some transmission of energy from the D7-brane 
into the D5-brane through the vertex. To set up such a configuration however one 
would need to carefully tune the transmission coefficients of each meson mode such 
that there was a stable groundstate (in the absence of a dynamical photon). 

Another interesting point to pursue further would be the addition of a magnetic 
field in the set up along the lines of [27, 41]. In this case one could investigate its 
effect on the chiral symmetry breaking and the behaviour of the spectral curves. We 
leave such investigations for future studies. 
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A. Alternative form of the Schrodinger potential 

Due to the numerical sensitivity in solving the eigenvalue problem we work with two 



- 30 - 



forms of the Schro dinger potential, one of which is most useful for calculating the 
mass spectrum and is described in the body of the text, while the second, shown 
here, is more useful for indicating qualitative features in the spectrum from the 
point of view of the potential. In the following parametrization we start again from a 
differential equation of the form (4.4) and transform it to a Schrodinger like equation 
along the lines of [42]. The transformed equation will have the following form 

dfy-v(()i> = o, (A.i) 

and the necessary steps that bring it into this form, starting from (4.4), are 

e c = p & y? = eif'^ip , (A. 2) 

with the potential given by 

V(() = -M^ + 1 -ff -l£ (A.3) 
Jo ^ Jo 4 Jo 

where 

f = e~<f, h = e^h. (A.4) 
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